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Abstract 

We study corrections to scaling in the 0(3)- and 0(4)-symmetric </> 4 
model on the three-dimensional simple cubic lattice with nearest neigh- 
bour interactions. For this purpose, we use Monte Carlo simulations in 
connection with a finite size scaling method. We find that there exists a 
finite value of the coupling A*, for both values of N, where leading cor- 
rections to scaling vanish. As a first application, we compute the critical 
exponents v = 0.710(2) and 77 = 0.0380(10) for N = 3 and v = 0.749(2) 
and T] = 0.0365(10) for N = 4. 



1 Introduction 



At a second order phase transition various quantities diverge with power laws. 
E.g. the magnetic susceptibility behaves as 

x~c ± \tr , a) 

where t = (T — T c )/T c is the reduced temperature. The subscripts + and — indi- 
cate the high and low temperature phase, respectively. 7 is the critical exponent 
of the magnetic susceptibility. 

The universality hypothesis says that for all systems within a given univer- 
sality class the exponent 7, as other critical exponents, takes exactly the same 
value. Note that the amplitudes C± depend on the details of the system, while 
the ratio C + /C_ is universal; i.e. it takes the same value for all systems within 
a universality class. A universality class is characterized by the spatial dimen- 
sion of the system, the range of the interaction and the symmetry of the order 
parameter. See ref. [[!]]. 

A central goal of the study of critical phenomena is to obtain precise estimates 
of universal quantities like the critical exponents or universal amplitude ratios. 

Often so called lattice spin models are used to study critical phenomena. The 
Ising model is the proto-type of such models. In three dimensions, where no 
exact solution of these models is available, the most precise results are obtained 
by the analysis of high temperature series or Monte Carlo simulations. For uni- 
versal quantities, a similar accuracy can be obtained with field-theoretic methods 
like the e-expansion or perturbation theory in three dimensions. For a detailed 
discussion see text-books on critical phenomena; e.g. refs. 0, ||. 

The accuracy of estimates of universal quantities extracted from high tem- 
perature series expansions or Monte Carlo simulations is limited by corrections 
to the power law ([[]): 

X = C ± |tr (1 + a± \t\ e + ...) , (2) 

where the correction exponent 9 = vuj pa 0.5 is universal. While the exponent 9 is 
universal, the amplitudes of corrections depend on the parameters of the system. 
Hence it is a very natural idea to search for parameters at which (non-analytic) 
corrections to scaling vanish. 

In the present work, we study the 0(iV)-invariant </> 4 model on simple cubic 
lattices with periodic boundary conditions. The lattices have the linear extension 
L in all three directions. The classical Hamiltonian is given by 

<xy> x x 

where <p x G M. N and < xy > denotes a pair of nearest neighbour sites on the 
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lattice. We study the canonical ensemble; the partition function is given by 



D[0] exp(-if) , (4) 

where j D [0] is a short hand for the N x L 3 dimensional integral over all com- 
ponents of the field at all lattice points. 

In the limit A = oo, the so called 0(iV)-invariant non-linear a model is recov- 
ered. In this limit the last term of the Hamiltonian forces the field to unit-length: 
cfl = 1. A = gives the exactly solvable Gaussian model. Along a critical line 
/3 C (A), the model undergoes a second order phase transition. For all A > 0, at 
given N, these transitions belong to the same universality class. Hence the criti- 
cal exponents, including the correction exponent 9, are the same for all values of 
A > 0. However, in eq. @ the correction amplitudes a± become functions of A. 
Hence, there is a chance to find a zero: a±(X*) = 0. Since leading corrections to 
scaling are due to a unique irrelevant scaling field w 3 , it follows that A* is unique 
for all quantities. 

Already in refs. [|], |5| it was suggested that for models that interpolate be- 
tween the Gaussian and the Ising model, there exists one value of the interpola- 
tion parameter, where leading corrections to scaling vanish. These models were 
studied with high temperature series expansions. 

More recently, this idea has been picked up in connection with Monte Carlo 
simulations and finite size scaling. In refs. [§, [F| a finite size scaling method has 
been proposed that allows to compute A* in a systematic way. While the details 
of the implementation of ref. || and |7j are different, the underlying physical 
idea is the same. 

This programme has been successfully implemented for N = 1 (Ising univer- 
sality class) H 0, H and N = 2 (XY universality class) |], [L0|. As a result, for 



these universality classes, the most precise estimates for critical exponents from 
Monte Carlo simulations were obtained. 

In a turn, the results for A* obtained from Monte Carlo have been used as 



input for the analysis of high temperature series of these models. In ref. [|TTJ the 
Ising universality class and in refs. |T0| , |I~2| the XY universality class was studied. 
As a result, the accuracy of the critical exponents was further improved and in 
addition, a large set of universal amplitude combinations has been obtained with 
unprecedented precision. 

In the present work we like to answer the question whether this programme 
can be extended to larger values of N. In ref. it was argued that for N > 3 
there exists no such A* for the Hamiltonian eq. (|3|) on simple cubic lattices. The 
argument is based on the large N expansion and the analysis of high temperature 
series expansions 
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The paper is organized as follows. In section 2 we briefly review the finite size 
scaling method that has been applied to eliminate corrections to scaling. The 
Monte Carlo algorithm is discussed in section 3. Details of the simulations are 
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given in section 4. In section 5 the numerical results are presented. Next we 
compare our results for the critical exponents with those given in the literature. 
Finally we give our conclusions. 



2 The finite size scaling method 
2.1 Phenomenological couplings 

Here, phenomenological couplings are used to locate A*. Nightingale, in his sem- 
inal paper on the phenomenological renormalization group [f[4"| |, proposed to use 
such a quantity to locate [3 C and to compute the renormalization group expo- 
nents. Phenomenological couplings are invariant under renormalization group 
transformations. Hence, they are well suited to detect corrections to scaling. 
The proto-type of such a quantity is the so called Binder-cumulant [fl~H]: 



U t = , (5) 

< m 2 > 2 



where 



m V 



h E & ( 6 ) 



is the magnetisation of the system. Note that also higher moments of the mag- 
netisation could be considered. 

Frequently the second moment correlation length divided by the linear ex- 
tension of the lattice ^nd/L has been studied. The second moment correlation 
length is defined by 

= \l A A Lt„lr\2 ' ( 7 ) 

where 

X = y ( [ i: <f>* ) ) (8) 




is the magnetic susceptibility and 




exp (i^j U J (9) 

is the Fourier-transform of the two-point correlation function at the lowest non- 
vanishing momentum. In order to reduce the statistical error, we averaged the 



results of all three directions of the lattice. Note that Nightingale 0] studied £/L, 
where £ is the exponential correlation length of a system of the size L D_1 x oo. 

The third quantity that we study is the ratio Z a /Z p of the partition function 
Z a with anti-periodic boundary conditions in one of the three directions and 
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Zp with periodic boundary conditions in all directions. Anti-periodic boundary 
conditions mean that the term Yl<xy> ' $v m ^ ne Hamiltonian is multiplied by 
— 1 for x = (Li,x 2 ,x 3 ) and y = (l,x 2 ,x 3 ). This ratio can be measured with the 
help of the boundary-flip algorithm, which is a version of the cluster algorithm. 



The boundary-flip algorithm was introduced in ref. |22[ for the Ising model. In 
[23|j the authors have generalized this method to the case of 0(iV)-invariant non- 
linear a models. As in ref. || we use a version of the algorithm that only measures 
Z a /Z p and does not perform the flip to anti-periodic boundary conditions. For a 



recent discussion of the algorithm see ref. |T0| 



2.2 The finite size method 

Below we shall briefly recall the theoretical basis of the finite size scaling method 
that has been developed in refs. || [7| §. See also ref. [|10 |. 



Let us denote the phenomenological couplings by R. We define a quantity R 
based on a pair of phenomenological couplings R\ and R 2 . 

First we define (3/ as the value of (3, where, at given A and L, R\ takes the 
fixed value Ri,f- 

R 1 (L,\,P f )=R 1J . (10) 

Hence f3f is a function of L and A. In the language of high energy physics, this 
is our "renormalization condition" . 
Next we define 

R(L,X)=R 2 (L,\,i3 f ) . (11) 

The behaviour of this quantity gives us direct access to corrections to scaling, as 
we shall see below. In the following we will frequently refer to R(L, A) as U R 2 at 

Let us consider linearized renormalization group transformations in the neigh- 
bourhood of the fixed point. In the following, we take explicitly into account the 
thermal scaling field u t and the leading irrelevant scaling field u%. In addition, 
we consider m 4 as representative of sub-leading irrelevant scaling fields. For the 
definition of scaling fields see e.g. the text-book [Q. 

In oder to derive the finite size scaling behaviour of R we block the lattice to 
a fixed size L'. The scaling fields u' t , u' 3 and u' A of the blocked system are given 
by 

u' t = a(X,P)L^ (12) 
with a(A, j3 c ) = by definition. The irrelevant scaling fields behave as 

u' 3 = uf\\)L- + 0((3-(3 c (X)) (13) 

and 

u' A = uf(\)L-+» + 0(P-p e {\)) . (14) 
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Here we have absorbed the constant factors V 1 ' v , L >ul and L /u}2 into a(X,j3), 
u^\X) and u\ (A), respectively. We always simulate at a vanishing external 
field. Hence also the corresponding scaling field Uh vanishes. While the value 
of ijj Rj 0.8 is well established (see e.g. [[UJ), the knowledge of uji is rather 
limited. Using the "scaling field method" , which is a version of Wilsons "exact 
renormalization group", the authors of ref. |T7]] find uj 2 = 1.78(11) for N = 3. 



Following the figure given in ref. [17| the value of u 2 for N = 4 is similar to 



that of iV = 3. In addition, there are corrections due to the breaking of the 



rotational invariance by the lattice. These corrections have an uo rot pa 2 [|18 
In the case of the Binder-cumulant and the second moment correlation length 
we have to expect corrections due to the analytic background of the magnetic 
susceptibility. This amounts to corrections with an exponent o>& a dfc = 2 — r] aa 2. 
In our numerical analysis, we will always pessimistically assume uj 2 = 1.6 when 
we estimate systematic errors due to sub-leading corrections. Since R keeps its 
value under renormalization group transformations we get 

R((3,X,L) = #««...) (15) 

Note that we are allowed to skip V from the list of arguments of R' since L' is 
fixed. Expanding R' in the scaling fields at the fixed point {u 1 - = 0} gives 

R = R* +D t u' t +D 3 u' 3 +D 4 u' A + ... (16) 
OR' 



with 



D; 



Inserting eqs. ([T^ , p!3|JT^) into eq. ([TB]) yields 



(17) 

K=o} 



R = R* + D t a(X,P f ) L 1/u + c(A) LT U + d(X) L^ 2 + ... , (18) 

where c(A) = D 3 uf \X) and d(X) = D^u^(X). Hence 

a(X, j3f) L l/v = D^ t (R 1>f - R{ - ci(A) L'" — di(X) L~ W2 ) + ... . (19) 

The subscript 1 has been introduced to indicate that the quantities are related 
to Ri. Taking R 2 at j3f gives 

R(L, X) = R* +c(X)L-" + d{X) L-" 2 + ... (20) 

with 

R* = R* + ^ (R hf - J?) , (21) 

c(A) = c 2 (A) - ^ Cl (A) (22) 
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and d is given by an analogous expression. 

In eq. fl20|) we have neglected corrections that are quadratic in the scaling 
fields. The most important of these is b c(A) 2 L~ 2uJ . The exponent 2u has 
roughly the same value as u^. However, since c(A) appears quadratically, we can 
safely ignore this term in the neighbourhood of A*. 

Eliminating leading corrections to scaling means in terms of eq. ( |20D to find 
the zero of c(A). One can imagine various numerical implementations to find this 
zero. Here we followed the strategy used in refs. || We have simulated the 
models close to the critical line for several values of A for various lattices sizes. 
The results are then fitted by an ansatz motivated by eq. (pCf) . The function c(A) 
is then approximated by interpolation between the A-values, where simulations 
have been performed. 

In the following, we will always use either Z a /Z p or ^nd/L as Ri and £7 4 as 
R 2 - Note that in ref. || we have only used Ri = Z a /Z p and in ref. || only 
Ri — i2nd/ L. Using both quantities gives us better control over systematic errors 
introduced by sub-leading corrections. 



3 The Monte Carlo Algorithm 

The 0(iV)-invariant non-linear a models were simulated with a cluster algorithm. 
For finite A, following Brower and Tamayo additional updates with a local 
Metropolis-like algorithm were performed to allow fluctuations of the modulus 
of the field <p x . Below, we give the details of the cluster algorithm and the local 
update. 



3.1 The wall-cluster update 

Several variants of the cluster algorithm have been proposed in the literature. The 
best know are the (original) Swendsen-Wang algorithm |20) and the single-cluster 
algorithm of Wolff f21fl . Here we have used the wall-cluster algorithm of ref. 0. 
In ref. J7J, for the 3D Ising model, a small gain in performance compared with 
the single-cluster algorithm was found. Such a comparison for N > 1 remains 
to be done. The main reason for using the wall-cluster algorithm here is that 
the wall-cluster update can be combined with the measurement of the ratio of 
partition functions Z a /Z p . In the wall-cluster update all clusters are flipped that 
intersect a plane of the lattice. 

The definition of a cluster is the same as in the Swendsen-Wang or the single- 
cluster algorithm. Following Wolff |2T|, only the sign of one component <px of 



the field is changed in a single step of the algorithm. This leads to the probability 

p rf = min[l,exp(-2/30i P Vl p) )] (23) 
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to delete the link < xy >. The links that are not deleted are frozen. A cluster is 
a set of sites that is connected by frozen links. 

In ref. [21] the component of the field parallel to a randomly chosen direction 
is used. Here we choose the 1 st , 2 nd , ... or N th component of <fi. This simplifies 
the implementation of the cluster update and the measurement of Z a /Z p . Also 
CPU-time is saved since in one update step only one component of the field 
has to be accessed. In order to compensate the algorithmic disadvantage of this 
restricted choice, we perform a global rotation of the field after a certain number 
of cluster updates. 



3.2 The local update of the 4 model 

We sweep through the lattice with a local Metropolis-type updating scheme. 
A proposal for a new field at the site x is generated by 

<j>f = 4$ + c ( r W - 0.5) , (24) 

where the are random numbers that are uniformly distributed in [0, 1), i runs 
from 1 to N. The proposal is accepted with the probability 

A = min[l, exp(-F' + H)} , (25) 

where H' is the Hamiltonian of the proposed field and H the Hamiltonian of the 
original field. The step-size c is adjusted such that the acceptance rate is about 
1/2. After this Metropolis step, we perform at the same site an over-relaxation 
step: 




Note that this step takes very little CPU-time. Hence it is likely that its benefit 
out-balances the CPU-cost. For lack of time, we did not carefully check this 
point. 



3.3 The update cycle 

Finally let us summarize the complete update cycle: 

• local update sweep 

(This step is omitted for the non-linear a models) 

• global rotation of the field 

• 3 x iV wall-cluster updates 

The sequence of the 3 x N wall-cluster updates is given by the wall in 1-2, 1- 
3 and 2-3 plane. In each of the three cases, an update is performed for all 
iV components of the field. For each of the three orientations of the wall we 
performed a measurement of Z a /Z p . 
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4 The simulations 



First we have simulated the 0(3)-invariant non-linear a model at the best esti- 
mate of p e = 0.693002(12) of ref. @. We used lattices of size L = 6, 8, 12, 
16, 24, 32 and 48. We have performed 25 x 10 6 measurements for L = 6 up to 
L = 32 and 10 7 measurements for L = 48. 

The N = 3 4 model was simulated at A = 2.0, 4.5 and 5.0 on lattices of the 
linear size L — 6, 8, 12, 16, 24 and 32. For A = 4.5 we simulated in addition 
L = 48. In all cases we performed 10 measurements. 

We have simulated the 0(4)-invariant non- linear a model at the estimate of 
(3 C = 0.935861(8) of ref. @. We studied lattices of size L = 6, 8, 12, 16, 24, 
32 and 48. We have performed 25 x 10 6 measurements for L = 6 up to L = 16, 
2 x 10 7 measurements for L = 24, 14 x 10 6 measurements for L = 32, and 105 x 10 5 
measurements for L = 48. 

The N = 4 4 model was simulated at A = 8.0, 12.0 and 14.0 on lattices of 
the linear size L = 6, 8, 12, 16, 24 and 32. For A = 12.0 we simulated in addition 
L = 48. As for N = 3, we performed 10 7 measurements for each parameter set. 

A measurement was performed after one update cycle. (See section |3~3| .) 

In all cases listed above, we performed 10 5 update cycles before measuring. 

In order to reduce the amount of data that is written to disc, we averaged 
during the simulation the results of 5000 measurements. These averages were 
saved. 

In the case of the 4 model we have simulated at estimates of /3f from Z a /Z p j, 
which were obtained from smaller lattice sizes that have been simulated before 
and/or short test-simulations. 

In our analysis the observables are needed as a function of (3. Given the 
large statistics we did not use the reweighting technique. Instead we used the 
Taylor-expansion up to the third order. The coefficients were obtained from the 
simulation. We always checked that the (3f are sufficiently close to the /3 of the 
simulation such that the error from the truncation of the Taylor-series is well 
below the statistical error. 

As random number generator we have used our own implementation of G05CAF 
of the NAG-library. The G05CAF is a linear congruential random number gen- 
erator with modulus m = 2 59 , multiplier a = 13 13 and increment c = 0. 

As a check of the correctness of the program and the quality of the random 
number generator we have implemented the following two non-trivial relations 
among observables: 

y.nn.x 
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and 

= ( ( E - ) ) - ^ E » (28) 

\ \j.nn.i / / y.nn.x 

where y.nn.x indicates that the sum runs over the six nearest neighbours of x. 
In order to enhance the statistics, we have summed eqs. fl27 , 28) over all sites 
x. We found that for all our simulations these equations are satisfied within the 
expected statistical errors. 

Most of the simulations were performed on 450 MHz Pentium III PCs. For 
our largest lattice size L = 48, one update cycle plus a measurement takes 0.67 
s, 0.58 s, 0.86 s and 0.75 s for the N = 3 4 model at A = 4.5, = 0.68622, 
the 0(3)-invariant non- linear a model at (3 = 0.693002, the N = 4 4 model at 
A = 12.0, (3 = 0.90843 and the 0(4)-invariant non-linear a model at /3 = 0.935861, 
respectively. 

In total, the whole study took about two years on a single 450 MHz Pentium 
III CPU. 



5 Analysis of the data 
5.1 Corrections to scaling 

In a first step of the analysis we have estimated ^nd/ L* and Z a /Z* by fitting the 
0(3)- and 0(4)-invariant non- linear a model data with the ansatz 

R(/3 C ) = R* + C L~ W , (29) 

where /3 C , R* and c are the parameters of the fit. We have fixed u — 0.8. As 
result we obtain Z a /Z* « 0.196 and ^nd/L w 0.564 for the 0(3) model and 
Z a /Z* w 0.1195 and £ 2n d/L* w 0.547 for the 0(4) model. In the following we 
shall use these numbers to set R±j. I.e. Z a /Z p j = 0.196 and ^nd/ 'Lf = 0.564 for 
N = 3 and Z a /Z pJ = 0.1195_and'«e 2n(i /L / = 0.547 for N = A. 

Next we have analyzed R to study corrections to scaling. To obtain a first 
impression, we have plotted our results for R with Z a /Z p j in figure |I]for N = 3 
and in figure |2| for N = 4. In both cases the range 6 < L < 32 is shown. 

Let us discuss in detail the N = 3 case. For the 0(3)-symmetric non-linear 
a model we clearly see an increase of R with increasing L over the whole range 
of lattices sizes. On the other side, for A = 2.0, R is decreasing. For A = 4.5 and 
A = 5.0 R stays almost constant. This behaviour suggests that leading corrections 
to scaling vanish at A* ~ 5. The behaviour of R for iV = 4 is qualitatively the 
same as for N = 3. The figure |2| indicates that A* ~ 13. 

In the following numerical analysis of the data we demonstrate that the be- 
haviour discussed above is indeed due to leading corrections to scaling and give 
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Figure 1: N = 3. The Binder-cumulant U at Z a /Z p j = 0.196 as a function of the 
lattice size L for A = 2.0, 4.5, 5.0 and oo. The dotted line should only guide the eye. 



an accurate estimate of A* and its error-bar. For this purpose we fitted our data 
for R with the ansatz 

R = R* + c(A) IT" (30) 

with R*, c(A) for each value of A and u as free parameters. For N = 3 and N = 4, 
we have performed such fits for three different sets of input data. These sets are 
given in table |I[ 

Our results for R* and u for N = 3 are given in table 0. 

The values of c from the same fits are summarized in table [5[ First of all 
we notice that our fit results for u are consistent with the estimates from field- 
theoretic methods. (See table [|) This indicates that the corrections are indeed 
dominated by the leading irrelevant scaling field U3. Given the statistical error 
and the variation of our result for u with the different data sets, we can not 
provide a more accurate estimate for u as the field-theoretic methods. 

Having convinced ourself that we really see leading corrections to scaling we 
extract an estimate for A* from the data given in table [3]. Therefore, we linearly 
extrapolate the results of c at A = 4.5 and A = 5.0. Taking the results from the 
set 1 we arrive at the estimate A* = 4.4(7), where the result from t7 4 at Z a /Z p j 
and U4 at £,2nd/Lf are consistent. The error-bar is computed from the variation 
of A* with the data-sets used for the fit. In the case of U4 at Z a /Z p j A* is roughly 
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Table 1: The sets of input data for fits with eq. For N = 3 as well as N = 4 
we have fitted with three different sets of input data. In column one the X-values 
for N = 3 and in column two the X-values for N = 4 are given. In columns three, 
four and five the lattice sizes L that have been included in the fits are listed. 



\, N = 3 


A, N = 4 


set 1 


set 2 


set 3 


2.0 


8.0 


16, 24, 32 


12, 16, 24 


8, 12, 16 


4.5 


12.0 


16, 24, 32, 48 


12, 16, 24, 32 


8, 12, 16, 24 


5.0 


14.0 


16, 24, 32 


12, 16, 24 


8, 12, 16 


oo 


oo 


16, 24, 32, 48 


12, 16, 24, 32 


8, 12, 16, 24 



Table 2: Results for R* and uj for N = 3 from fits with the ansatz ft30{). The 
data-sets that have been used for the fits are given in table |J. The results for c 
from the same fits are summarized in table I 



set 


X 2 /d.o.l 


R* 




U± ai 


' Z a /Z p j 


1 
2 
3 


2.32 
1.79 
1.14 


1.14018(9) 
1.14028(9) 
1.14021(8) 


0.743(23) 
0.749(18) 
0.799(13) 


U 4 ai 


C (,2nd / L 


1 
2 
3 


2.16 
0.90 
4.06 


1.13977(10) 
1.14000(11) 
1.14038(10) 


0.741(22) 
0.732(17) 
0.775(12) 



Table 3: Results for c for N 
R* and uj are given in table [^. 



3 from fits with the ansatz (3i). The results for 



set 


5(2.0) 


c(4.5) 


c(5.0) 


c(oo) 


U4 at Z a /Zpj 


1 


0.0392(33) 


-0.0008(10) 


-0.0042(8) 


-0.0390(24) 


2 


0.0394(25) 


-0.0011(8) 


-0.0049(7) 


-0.0408(17) 


3 


0.0451(17) 


-0.0007(6) 


-0.0053(5) 


-0.0461(11) 


U4 at (,2nd/ Lf 


1 


0.0464(37) 


-0.0007(11) 


-0.0056(9) 


-0.0464(26) 


2 


0.0436(27) 


-0.0026(8) 


-0.0074(7) 


-0.0475(18) 


3 


0.0452(18) 


-0.0064(6) 


-0.0115(5) 


-0.0566(12) 
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Figure 2: N = 4. The Binder-cumulant {/ at Z a /Z p j = 0.1195 as a function of the 
lattice size L for A = 8.0, 12.0, 14.0 and oo. The dotted line should only guide the eye. 



the same for all three sets. However for U4 at ^2nd/Lf there is a drift to larger 
results for A* as the lattice sizes increase. Assuming a convergence proportional 
£-w 2 +w ~ £-°- 8 we arrive at our error-estimate. 

In tables [| and |] we have summarized our results for R* and u and c for 
N = 4. As for JV = 3 we see that the values for u are consistent with the 
field-theoretic results. 

Clearly, the sign of c(oo) is negative and that of c(8.0) is positive. Hence there 
exists a A* with c(A*) = 0. From linear interpolation of the result for A = 12.0 
and A = 14.0 we arrive at A* = 12.5(4.0). The error-bar has been computed in 
the same way as for N = 3. 

We have also used a slightly different approach to compute uj from the avail- 
able data. We have analyzed the difference of R at A = 2.0 and A = 00. Obviously, 
R* is cancelled this way. In addition one might expect that corrections that are 
quadratic in the scaling field M3 cancel since |c(2.0)| and |c(oo)| are almost the 
same. In addition, the analysis of numerical data for iV = 1 in ref. || and N = 2 
in ref. suggests that sub-leading corrections also cancel to a large extent. 
Therefore, we have fitted our data with the ansatz 

R(L,\)\ x=20 -R(L,\)\ x=oo = AcL-" . (31) 

Our results for N = 3 and the corresponding results for iV = 4 are given in table 
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Table 4: Results for R* and uj for N = 4 from fits with the ansatz (3C). The 
data-sets that have been used for the fits are given in table 0. The results for c 
from the same fits are summarized in table [|. 



set 



xVd.o.f. 



R* 



to 



U 4 at Z a /Z p , 



f 



1 

2 
3 



0.75 
0.84 
2.07 



1.09441(6) 
1.09450(7) 
1.09466(6) 



0.798(52) 
0.748(41) 
0.761(27) 



U± at ^2nd/ L 



1 

2 
3 



0.77 
1.08 
7.51 



1.09458(7) 
1.09474(8) 
1.09522(8) 



0.761(50) 
0.735(37) 
0.761(25) 



Table 5: Results for c for N - 
R* and uj a re given in table [i . 



4 from fits with the ansatz (\3Q). The results for 



set 


c(8.0) 


c(12.0) 


c(14.0) 


c(oo) 


U4 at Z a /Z p j 


1 


0.0091(19) 


0.0010(8) 


-0.0023(7) 


-0.0160(22) 


2 


0.0073(12) 


-0.0001(6) 


-0.0021(5) 


-0.0145(14) 


3 


0.0061(7) 


-0.0013(4) 


-0.0033(3) 


-0.0165(8) 


U4 at £,2nd/ Lf 


1 


0.0083(18) 


0.0011(8) 


-0.0035(6) 


-0.0183(23) 


2 


0.0065(12) 


-0.0017(6) 


-0.0040(5) 


-0.0184(15) 


3 


0.0030(7) 


-0.0059(4) 


-0.0080(3) 


-0.0239(10) 
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Table 6: Results for 00 for N = 3 and N = 4 obtained by fitting with the 
ansatz (I31U. In the first column we give N and the phenomenological coupling 
that has been used to determine (3f. In the fit all lattice size with L min < L < 
L mnr have been taken into account. 





Lmin 


Lmax 


X 2 /d.o.f. 


us 


N 


= 3, U A at Z a /Z p = 0.196 


6 


32 


1.67 


0.796(7) 






8 


32 


0.73 


0.781(10) 


N 


= 3, U A at £ 2 /L = 0.564 


6 


32 


0.65 


0.769(6) 






8 


32 


0.65 


0.766(9) 


N 


= 4, U A at Z a /Z p = 0.1195 


6 


32 


0.54 


0.780(15) 






8 


32 


0.40 


0.765(22) 


N 


= 4, U 4 at C2/L = 0.547 


6 


32 


0.40 


0.774(14) 






8 


32 


0.38 


0.764(20) 



|^. In the case of iV = 4 we have taken the difference of R at A = 8.0 and A = 00. 

First we notice that a x 2 /d.o.f. ~ 1 is already reached for L min = 6. For 
such a small L min , fits with ansatz (|30|) produce x 2 /d.o.f. = 2.8 for N = 3, 
at Z a /Z p = 0.196 and xVd.o.f. = 13.3 for iV = 3, U A at f 2 /L = 0.564. This fact 
indicates that the above mentioned cancellations indeed occur. 

The results for u are consistent with those of the field theoretic methods. 
Certainly our approach is very promising to give competitive results for the cor- 
rection exponent u. However, we would like to have still larger statistics and 
a larger range of lattice sizes to give a sensible estimate of the systematic error 
caused by sub-leading corrections. 



5.2 Critical exponents 

We have computed the critical exponents u and n using well established finite 
size scaling methods. Below we shall only discuss in detail the results for N = 3. 
The analysis for N = 4 has been performed analogously. 

The exponent v is computed from the slope of a phenomenological coupling 
R at p f (see eq. (0)): 

9 Ar aL1 "- (32) 

As it was pointed out in |2H}| , replacing (3 C by (3f simplifies the error-analysis, 
since the error in j3 c needs not to be propagated. 

In our study we have considered three different choices of R. Hence in eq. (j32j) 
we could in principle consider 9 different combinations. Below we shall restrict 
ourself to six choices: j3f is fixed either by Z a /Z p j or fand/Lf- AVe consider the 
slope of all three phenomenological couplings R. 
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Table 7: Results for v from fits with the ansatz (p^j. The range of lattice sizes 
is always L m i n = 16 and L max = 32. We have used six different combinations of 
R. These are given by CI: The slope ofU^ at Z a /Z p j; C2: The slope of Z a /Z p 
at Z a /Z p j; C3: The slope of^nd/L at Z a /Z p j; C4: The slope ofU^ at £,2nd/Lf', 
C5: The slope of Z a /Z p at £ 2n d/Lf; C6: The slope of&nd/L at &nd/L f . 



A 


CI 


C2 


C3 


C4 


C5 


C6 


2.0 


0.7164(23) 


0.7088(10) 


0.7136(13) 


0.7181(24) 


0.7106(12) 


0.7135(12) 


4.5 


0.7071(21) 


0.7083(9) 


0.7115(12) 


0.7074(21) 


0.7086(10) 


0.7115(12) 


5.0 


0.7078(21) 


0.7085(9) 


0.7114(11) 


0.7071(21) 


0.7077(10) 


0.7114(11) 


oo 


0.7076(12) 


0.7127(5) 


0.7142(6) 


0.7056(12) 


0.7136(13) 


0.7141(6) 



Table 8: Results for v from fits with ansatz ( |3^ ) for N = 3 at A = 4.5; always 



Lmax = 48. The combinations C1,...,C6 are explained in the caption of table |/f. 



Lmin 


CI 


C2 


C3 


C4 


C5 


C6 


8 

12 
16 


0.7113(7) 

0.7111(10) 

0.7101(13) 


0.7087(3) 
0.7100(4) 
0.7099(6) 


0.7150(4) 
0.7132(5) 
0.7120(7) 


0.7102(7) 

0.7107(10) 

0.7102(14) 


0.7075(3) 
0.7096(5) 
0.7100(7) 


0.7150(4) 
0.7132(5) 
0.7120(7) 



First we like to study how much the result of v from fits with eq. ( |32"D depends 
on leading corrections to scaling. For this purpose we have fitted for L min = 16 
and L max = 32 our data for all available values of A for all six combinations of 
R. The results of these fits are summarized in table 0. 

The variation of the results with A are rather small. The largest variation 
we find for the combination CI, where we get v = 0.7164(23) for A = 2.0 and 
v = 0.7076(12) for A = oo. Hence, we expect that for A = 4.5 the effect of 
corrections to scaling should be smaller than 0.0001. 

Next, let us discuss in more detail the results for A = 4.5 which is closest to 
our estimate of A*. Results from fits with ansatz (|32|) are summarized in table [8[ 
We see that the results approach each other as L min is increased. In the case of 
the slope of U4, the results for v stay almost constant as L min is varied. For the 
slope of Z a /Z p we see a slight increase of the estimate of v. On the other hand 
for the slope of ^nd/L we see a decrease. Assuming that this behaviour is caused 
by the sub-leading corrections, we conclude that v = 0.7120 from ^nd/L is an 
upper bound. Given the larger stability of the estimate from Z a /Z p we take as 
our final estimate v = 0.710(2). 

Next we computed the exponent n from the finite size behaviour of the mag- 
netic susceptibility 

XL = cL^ . (33) 
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We restrict the discussion to Z a /Z p j since fixing j3f by ^nd/Lf gives very similar 
numbers. 

For the estimate of rj we see a much stronger dependence on A than for v. 
Fitting with ansatz (|33"D and L min = 16 , L max = 32 we obtain r] = 0.03982(33), 
0.03614(30), 0.03541(31) and 0.03269(18) for A = 2.0, 4.5, 5.0 and oo, respec- 
tively. 

Fitting the data at A = 4.5 with L max = 48 yields r) = 0.03581(14), 0.03668(19) 
and 0.03736(32) for L min = 12, 16 and 24, respectively. For L min = 12 we get 
X 2 /d.o.f.= 18.5. The strong dependence of the result on L min and the large 
X 2 /d.o.f. at L min = 12 indicates that there are sizeable sub-leading corrections. 

Fitting with an ansatz that includes an analytic background term 

X \p f = c L 2 ~" + b (34) 

yields x 2 /d-o.f.= 0.38 already for L min = 8. The results are r] = 0.0384(2), 
0.0386(4) and 0.0381(6) for L min = 8, 12 and 16, respectively. To see the effect 
of leading corrections on this fit we have in addition fitted the data for A = 5.0 
with L min = 8 and L max = 32. We get rj = 0.0378(4). 

As our final estimate we quote rj = 0.0380(10). The error-bar takes into 
account statistical errors as well a systematic errors due to leading and sub- 
leading corrections. These errors are estimated from the spread of the results of 
the various fits discussed above. 

With a similar analysis we arrive at v — 0.749(2) and rj = 0.0365(10) for the 
0(4) universality class. 



6 Comparison with results from the literature 



Here we like to compare our results for the critical exponents v and rj with 
previous Monte Carlo studies of the 0(3)- and 0(4)-invariant non-linear a models. 
In addition we give selected results from high temperature series, perturbation 
theory in three dimensions, and the e-expansion. The results are summarized in 
tables |^ and [It] for N = 3 and N = 4, respectively. For more references on field 
theoretic methods and other methods, not discussed here, see e.g. ref. fl~6|] . 

All Monte Carlo studies listed in the tables § and [H| use a simple cubic lattice. 
In addition, in ref. [29| the body centred cubic lattice is studied. (In table |we 
only give the sc results.) In all studies the lattice sizes are smaller or equal 
L = 48, except for ref. where in addition L = 64 is simulated. The results of 
the MC studies cited above are extracted from ansatze like eqs. (|32]j33| ). (Mostly 
(5 C is used instead of j3f.) We see that almost all Monte Carlo results for v are 
consistent with ours. On the other hand, the results for rj are systematically too 
small, except for ref. [2~1|. 

This behaviour can be well understood with our results of section |5.2| . The 
estimates for rj from the ansatz (|3"B|) are clearly effected by corrections to scaling. 
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Our results from the 0(3)- and 0(4)-invariant non-linear a models for rj are 
systematically lower than our final results from the 4 models at A*. On the 
other hand, the results for v given in table [7| show only little variation with A; 
i.e. little dependence on leading corrections to scaling. 



The authors of ref. [p^|, who for the first time tried to take into account 
leading corrections to scaling in the analysis of their data, arrive at rather similar 
conclusions how leading corrections to scaling affect the estimates of rj and v. 
They extrapolated their results for rj assuming L~ w corrections. 



However from our analysis of section |572| we know that the estimates of 77 
obtained from lattices with L < 48 are also strongly affected by sub-leading 
corrections with LO2 ~ 2. 

Hence, extrapolating only in L~ u leads to a wrong amplitude for the L~ w 



corrections. As result, the final estimate of ref. |24| for 77 is too large compared 
with our result. 

There exists a large number of refs. on the e-expansion and the perturbative 
expansion in three dimensions in the literature. As an example, we have chosen 
the result of a recent analysis by Guida and Zinn- Justin ||16|| . We notice that 
these results are consistent with ours for rj and v within the quoted errors. Also 
note that the error-bars of our estimates for v and rj are smaller than those of 
the field-theoretic estimates. 

There exists also a number of publications on the analysis of high temperature 
series. In the tables we give the results of a recent analysis [0 using inhomoge- 
neous differential approximants. The coefficients of the high temperature series 
of x an d /i2 are computed up to (3 21 . In our tables, we only give the results from 
the unbiased analysis of the simple cubic lattice series. In addition the authors 
analyze the body centred cubic lattice. The authors also give results obtained 
from a so called ^-biased analysis, where they make use of the numerical results 
for 9 = ujv obtained from field-theoretic methods. It is interesting to notice that 
these baised results (not given in our tables) tend to be less consistent with our 
results than the un-biased results which we quote in tables and |TU[ 

In tables ^| and 10|, we give for uj the average of the result from Z a /Z p j and 



£,2nd/Lf with L min = 8 taken from table |6|. We make no attempt to estimate 
the systematic errors due to sub-leading corrections. Certainly these errors are 
larger than that quoted for the field-theoretic estimates of ref. ||16|| . It is however 



interesting to note that our results are consistent with those of ref. [16 



Also the Monte Carlo result of ref. [241 for N = 3 is consistent with ours. 



However we cannot confirm their surprising result for N = A. 

7 Conclusions 

In this study we have demonstrated that the programme of refs. 0, [7| to eliminate 
leading corrections to scaling in the three dimensional </> 4 model can be extended 
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Table 9: Results for the critical exponents of the 0(3) universality class from 
various methods. The numbers for rj that are marked by a * are computed from 
7/V = 2 — rj. Details are discussed in the text. 



method 


ref. 


V 


V 


UJ 


IMC 


present work 


0.710(2) 


0.0380(10) 


0.773 


MC 




m 


0.7128(14) 


0.0413(15)(1) 


0.78(2) 


MC 




m 


0.642(2) 


0.020(1) 




MC 




m 


0.704(6) 


0.027(2) 




MC 




m 


0.704(6) 


0.028(2) 




MC 




pi 


0.706(9) 


0.031(7) 




MC 




f^9 


0.7036(23) 


0.0250(35) 




HT 


N 


0.715(3) 


0.036(10)* 




d = 3 PT 


[16 


0.7073(35) 


0.0355(25) 


0.782(13) 


e-expansion 




16 


0.7045(55) 


0.0375(45) 


0.794(18) 



Table 10: Results for the critical exponents of the 0(4) universality class from 
various methods. The numbers for rj that are marked by a * are computed from 
7/V = 2 — rj. Details are discussed in the text. 



method 


ref. 


V 


V 


UJ 


IMC i 


jresent worl 


c 0.749(2) 


0.0365(10) 


0.765 


MC 


|24| 


0.7525(10) 


0.0384(12) 


1.8(2) 


MC 


|31| 


0.7479(90) 


0.0254(38) 




HT 


[30| 


0.750(3) 


0.035(9)* 




d = 3 PT 


|16| 


0.741(6) 


0.0350(45) 


0.774(20) 


e-expansion 


m 


0.737(8) 


0.0360(40) 


0.795(30) 
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to N = 3 and N = A. In particular, we have found A* = 4.4(7) for N = 3 
and A* = 12.5(4.0) for N = 4. Based on this result we have computed the 
critical exponents v and 77 from finite size scaling. In particular in the case of rj, 
the error-bar could be reduced considerably compared with previous Monte Carlo 
simulations or field theoretic methods and the analysis of high temperature series. 

Since the CPU-time that was used for the present study is still moderate, 
further progress can be made by just enlarging the statistics and simulating larger 
lattices. 

Also, our results for A* can be used as input for the analysis of high temper- 



ature series analogous to refs. [11, 12 



The principle question raised in ref . jTTJ , whether the programme to eliminate 



leading corrections is restricted to N < N c , where N c is finite, remains open. 
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